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Abstract 



We have developed an image restoration technique based on the Richardson- Lucy algorithm optimized 
for GLAST-LAT image analysis. Our algorithm is original since it utilizes the PSF (point spread 
function) that is calculated for each event. This is critical for EGRET and GLAST-LAT image 
analysis since the PSF depends on the energy and angle of incident gamma-rays and varies by more 
than one order of magnitude. EGRET and GLAST-LAT image analysis also faces Poisson noise 
due to low photon statistics. Our technique incorporates wavelet filtering to minimize noise effects. 
We present studies of EGRET sources using this novel image restoration technique for possible 
identification of extended gamma-ray sources. 
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1. INTRODUCTION 



The EGRET (Energetic Gamma Ray Experiment Telescope) instrument onboard the CGRO (Compton Gamma 
Ray Observatory) detected 271 gamma-ray sources in the energy band 0.1-10 GeV.^ However, 170 of these 
sources do not have clear counterparts at other wavelength and many of them still remain unidentified. Most 
of unidentified EGRET sources are located at low Galactic latitude and are considered to be of Galactic origin. 
Although some of them can be extended sources such as supernova remnants, it is difficult to identify extended 
sources due both to poor PSF (point-spread function) of the EGRET and to the presence of the Galactic diffuse 
gamma-ray background. An earlier attempt to improve the EGRET image using Richardson-Lucy deconvolution 
technique^ with wavelet denoisingby Charalabides et al.^ suffered from a rapidly varying and energy dependent 
PSF. We have developed a novel image restoration technique that incorporates the PSF calculated for each event 
as described in the following section. This technique will be critical for analysis of extended sources in the LAT 
(Large Area Telescope) instrument of the GLAST (Gamma-ray Large Area Telescope) since its PSF varies by 
two order of magnitude over the wide energy band, 0.1-300 GeV. 



2. IMAGE RESTORATION TECHNIQUE 

The observed image (j){x) is a convolution of the true image '0(^) and the instrument response function P(x|^), 
where P{x\^) is the probability that the photon is observed at x when the true position is ^. Based on Bayesian 
approach, the original image '0(0 can be obtained iteratively: 
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where x^ is the observed position of the kth event and P/c(^/c|0 probability to observe the event at x^ 

when the true position is ^. (The first formula is generalized for the case where the PSF is assigned for each 
event, while the second formula in the square brackets is the original version shown as a reference.) It is proven 
mathematically that '^^^^(O yields larger likelihood than at each iteration. However, it suffers from the 

amplification of the noise (mostly Poisson noise due to poor statistics in the EGRET analysis) as the number of 
iterations is increased. Efforts have been made to define the objective criteria to stop the iteration before such 
effects become prominent.^ Wavelet denoising can be used to suppress the noise in the residual between the 
observed and expected image in each iteration.^ The residual, p^(x), is defined as p^(x) = 0(x) — 0^(x), where 
0^(x) is the convolved image defined as: 
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Here n(x|C) is the weighted average of the instrument response at a given position x. The residual is decomposed 
to J wavelet scales as p^{x) = cj{x) + X]j=i ^j(^)7 where cj{x) is the last smoothed image and w denotes each 
wavelet scale. The residual can be filtered based on its significance as: 



p'{x) = cj{x) + ^Mj{x)wj{x), 



(3) 



where Mj{x) is if the Wj{x) is consistent with noise and 1 otherwise. We currently use the 99% Poisson 
probability as the consistency threshold. The filtered residual is incorporated in the iteration process using the 
filtered image, (x) = p'^{x) + (j/{x): 
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The generalized formulas described here are completely equivalent to the non-generalized versions if Pk{x\£) is 
common for all events. 
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When we have known point sources, the probabihty density function, can be decomposed into two 

components: '^^(0 = i^psiO + V^dfIO? where V^p^(^) represents a point source component (mostly zero apart 
from (5-functions at the positions of the point sources) while '^^^(O represents a diffuse component.^ Wavelet 
filtering is applied only to the diffuse component. This dual channel method enables us to incorporate sharp 
point sources while keeping the smooth diffuse component. This effectively works as point source subtraction 
without negative values in '^^^(O (straight point source subtraction often causes negative values due to noise). 



3. APPLICATION TO EGRET IMAGE ANALYSIS 

In this section, we demonstrate our image restoration technique by analyzing three fields observed by the EGRET. 
The PSF is calculated for each event as a convolution of PSFs in different energy bins, according to the energy 
dispersion of the given energy measurement. Only photons with energy greater than 0.5 GeV are used in this 
analysis in order to minimize the processing time without sacrificing the image quality. Lower energy photons do 
not contribute to improving the quality of the restored image. (It does not degrade the image quality either since 
such photons are weighted appropriately via large PSF). The results presented in this paper are still preliminary 
and should not be used for any quantitative analysis. 

3.1. Large Magellanic Cloud 

The LMC (Large Magellanic Cloud) is one of few extended sources resolved by EGRET. Figure 1 shows a 
smoothed photon count map of the LMC in the left panel, a restored image of the LMC in the middle panel, 
and the IRAS LMC image overlaid with the contour of the restored EGRET image in the right panel. The 
restored EGRET image is apparently sharper with better contrast than the original count map while avoiding 
enhancement of spurious structures (sometimes suppressing them) in the background (The color scale is common 
for both images). The bright spot of the restored image is well correlated with a dense region indicated by the 
IRAS image, which traces dusty interstellar clouds and associated star-forming regions. 
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Figure 1. Smoothed count map (left), restored image (middle) of the LMC by EGRET, and the IRAS LMC image overlaid with 
contour of the restored EGRET image (right). 



3.2. 3EG J1234-1318 

We also analyzed several unidentified EGRET sources at high Galactic latitude where the background diffuse 
component is very low. We found one source that is inconsistent with one point source as shown in Figure 2. The 
left panel shows the smoothed count map of the 3EG J1234-1318 while the right panel shows the restored image. 
A white cross and a contour in each image indicate the source location and its 90% confidence level contour in 
the original 3EG catalog (The 3EG catalog is produced using photons with £^ > 0.1 GeV). The restored image 
has much better contrast and clearly shows elongated structure indicating the possibility of two point sources 
with a distance of ~1.5°or an extended source. Since the Galactic latitude of this source is ~50° it is very likely 
to be two extra- Galactic point sources. If it is the case, new point source position would be preliminarily (RA, 
DEC)=(187.9°, -14.1°) for the brighter one and (188.9°, -13.0°) for the dimmer one instead of the original 
position, (188.5°, —13.3°). Another possibility is that this is due to a nearby molecular cloud interacting with 
cosmic-rays, which could be extremely useful to obtain cosmic-ray fiux (spectrum and intensity) in that region. 
This result demonstrates that our image restoration technique can be quite effective in identification of extended 
sources and source confusions. 
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Figure 2. Smoothed count map (left) and restored image (right) of the 3EG J1234-1318 by EGRET. 

3.3. Galactic Plane 

Although the Milky Way can be resolved by the EGRET due to its large size, smaller scale structure in the 
Galactic plane is not easy to resolve due to large EGRET PSF and contamination of the point sources. Our 
image restoration technique can reduce these effects as shown in Figure 3. The left panel shows the smoothed 
count map, the middle panel shows the restored image with known point sources included by the dual channel 
method, and the right panel shows the restored image of only the diffuse channel from the previous image, which 
is equivalent to the point source removal. The final EGRET image of the Galactic plane has smooth background 
and less structures associated with "point sources" that may or may not be real point sources. 




Figure 3. Smoothed count map (left), restored image with known point sources included by the dual channel method (middle), 
and restored image of the only diffuse channel (right) of the Galactic plane around the Galactic center by EGRET. 

4. CONCLUSIONS 

Generalization of a widely-used image restoration technique is introduced to incorporate the rapidly varying PSF 
in the EGRET and GLAST-LAT although it is not trivial since the wavelet filtering is performed on pixelated 
images. We demonstrated the effectiveness of our technique using three fields observed by the EGRET. Our 
technique is confirmed to be useful to obtain sharper images of extended sources, to identify source confusions, 
and to cleanly remove point sources from the image of extended sources. 
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